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Abstract. We examine the properties of Wolf-Rayet (WR) stars predicted by models of rotating stars taking 
account of the new mass loss rates for O-type stars and WR stars (Vink et al. l20UOII2U01l Nugis & Lamers [20001 
and of the wind anisotropies induced by rotation. We find that the rotation velocities v of WR stars are modest, 
i.e. about 50 km s~^, not very dependant on the initial v and masses. For the most massive stars, the evolution of 
V is very strongly influenced by the values of the mass loss rates; below ^^12 Mq the evolution of rotation during 
the MS phase and later phases is dominated by the internal coupling. Massive stars with extreme rotation may 
skip the LBV phase. 

Models having a typical v for the O-type stars have WR lifetimes on the average two times longer than for 
non-rotating models. The increase of the WR lifetimes is mainly due to that of the H-rich eWNL phase. Rotation 
allows a transition WN/WC phase to be present for initial masses lower than 60 Mq. The durations of the other 
WR subphases are less affected by rotation. The mass threshold for forming WR stars is lowered from 37 to 22 
Mq for typical rotation. The comparisons of the predicted number ratios WR/O, WN/WC and of the number 
of transition WN/WC stars show very good agreement with models with rotation, while this is not the case for 
models with the present-day mass loss rates and no rotation. As to the chemical abundances in WR stars, rotation 
brings only very small changes for WN stars, since they have equilibrium CNO values. However, WC stars with 
rotation have on average lower C/He and O/He ratios. The luminosity distribution of WC stars is also influenced 
by rotation. 
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$H ■ 1- Introduction loss rates, is no longer applicable at present. Indeed, the 

, estimates of the mass loss rates during the WR phase are 

Wolf-Rayet stars play an important role in Astrophysics, reduced by a factor of 2 to 3, when account is taken of 

as signatures of star formation in galaxies and starbursts, ^^le clumping effects in the wind (Nugis and LamersEHnni 

as injectors of chemical elements and of the radioactive iso- Hamann and Koesterke [T!?99.1 . Also, the mass loss rates 
tope 26 Al, as sources of kinetic energy into the interstellar Q-type stars have been substantially revised (and in 

medium and as progenitors of supernovae and maybe as general reduced) by the new results of Vink et al. 



progenitors of 7-ray bursts. As rotation affects all outputs ^^o also account for the occurrence of bi-stability 
of stellar models (Meynet & Maeder 12009 , the main pur- ^^^^^^^ ^^^^^ ^^^^^^ ^j^g ^-^^^ properties and mass loss 
pose of the present paper is to examine some consequences ^i^^g context, it is quite clear that with these 
of rotation on the properties of Wolf-Rayet (WR) stars. ^^ass loss rates the predicted numbers of WR stars 
Let us recall some difficulties faced by non-rotating by standard non-rotating models would be much too low 
stellar models of WR stars. A good agreement between with respect to the observations, 
the predictions of the stellar models for the WR/O num- 
ber ratios and the observed ones at different metallicities A second difficulty of the standard models with mass 
in regions of constant star formation was achieved pro- loss was the observed number of transition WN/WC stars, 
vided the mass loss rates were enhanced by about a factor These stars simultaneously show some nitrogen character- 
of two during the MS and WNL phases (Meynet & al. istics of WN stars and some carbon of the further WC 
I1994|l . This solution, which at that time appeared reason- stage. The observed frequency of WN/WC stars among 
able in view of the uncertainties pertaining to the mass WR stars is around 4.4 % (van der Hucht i2001fl . while 
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the frequency predicted by the standard models without 
extra mixing processes are lower by 1-2 orders of magni- 
tude (Maeder & Meynet A third difficulty of the 
standard models as far as WR stars were concerned was 
that there were relatively too many WC stars with respect 
to WN stars predicted (see the review by Massey [5D03). 
These difficulties are the signs that some process is missing 
in these models. 

We explore here the effects of rotation on the WR stars 
and their formation. Rotation already improved many pre- 
dictions of the stellar models, bringing better agreement 
with the observations. In particular, star models with ro- 
tation have at last enabled us (Maeder & Meynet ESnj 
to solve the long-standing problem of the large number 
of observed red supergiants at low metallicities Z in the 
SMC, while current models with any kind of reasonable 
mass loss usually predicted almost no red supergiants at 
low Z. Also, in addition to several other points, the new 
models also show significant He- and N-enhancements in 
0-type stars on the MS (Heger and Langer [?UUUI Meynet 
& Maeder r2000(l . as well as in supergiants at various metal- 
licities, as required by observations. 

Sect. 2 summarizes the physics of these models. The 
evolution of the inner and surface rotation is discussed in 
Sect. 3. The evolutionary tracks are shown in Sect. 4. The 
effects of rotation on the WR star formation process and 
on the WR lifetimes are discussed in Sect. 5. Comparisons 
with the observed WR populations are performed in Sect. 
6. Finally the predicted abundances in WR stars are dis- 
cussed in Sect. 7. 

2. Physics of the models 

The present grid of models at solar metallicity is based in 
general on the same physical assumptions as in the grid by 
Meynet & Maeder CUDD", paper V). However, in addition 
we include here several improvements that have appeared 
since. These are: 

— As reference mass loss rates in the case of no rotation, 
we use the recent data by Vink et al. 1^11117111^71)1 . A 
bi-stability limit intervenes at Tcff = 25000 K and pro- 
duces some steep changes of the mass loss rates during 
MS evolution. For the domain not covered by these au- 
thors we use the results by dc Jager et al. H1988|l . Since 
the empirical values for the mass loss rates employed 
for non-rotating stars are based on stars covering the 
whole range of rotational velocities, we must apply a 
reduction factor to the empirical rates to make them 
correspond to the non-rotating case. The same reduc- 
tion factor as in Maeder & Meynet l|2001|l was used 
here. 

— During the Wolf-Rayet phase we apply the mass loss 
rates proposed by Nugis and Lamers H2000fl . These 
mass loss rates, which account for the clumping ef- 
fects in the winds, are smaller by a factor 2-3 than 
the mass loss rates used in our previous non-rotating 
stellar grids (Meynet et al. ll994|l . 



— The significant increases of the mass loss rates M 
with rotation are taken into account as explained in 
Maeder & Meynet (2000a). This treatment, in contrast 
to the one usually used (Friend and Abbott [l98f)|l . ac- 
counts for the non-uniform brightness of rotating stars 
due to the von Zeipel theorem and for the fact that 
the Eddington factor is also a function of the rota- 
tion velocity. Due to the non-uniform brightness of 
a rotating star, the stellar winds are anisotropic. The 
rotation-induced anisotropy of the stellar winds is also 
accounted for during the Main Sequence (MS) phase 
according to the prescription of Maeder H2002|l . For 
0-type stars, the polar winds are enhanced and this 
means that these stars lose a lot of mass, but a mod- 
erate amount of angular momentum. For stars cooler 
than Tcff — 25000 K, the equatorial ejection becomes 
significant and leads to a larger ejection of angular 
momentum. After the MS phase, the ratios of the sur- 
face angular velocity to the break-up velocity are in 
general much too low for the anisotropics to be impor- 
tant. Interestingly enough, in cases of extreme mass 
loss such as in WR stars, the timescale for the trans- 
mission of the extra torque applied at the stellar sur- 
face by the anisotropics is longer than the timescale 
for the mass removal of the considered layers, so that 
the effects of the anisotropics of the stellar winds on 
the internal rotation are rather limited f Maeder ^ Ml^ . 

— In paper V, we did not account for the effects of hori- 
zontal turbulence on the shears (Talon & Zahn [TM7)l . 
We account for it in the present work in the same man- 
ner as we did in paper VII (Maeder & Meynet I2()()l|l . 
Let us recall that the horizontal turbulence, expressed 
by a diffusion coefficient Dh , tends to reduce the shear 
mixing in regions of steep //-gradients, making the 
diffusion coefficient I?shcar proportional to Z?h rather 
than to the thermal diffusivity. The difference can 
be large and it leads to smaller surface enhancements 
in the products of CNO burning. In regions of low 
/x-gradients, the horizontal turbulence moderately in- 
creases the shear diffusion coefficient, making it pro- 
portional to -Dh + K rather than to K only, with little 
consequence for the surface enrichments. 

— A moderate overshooting is included in the present ro- 
tating and non-rotating models. The radii of the con- 
vective cores are increased with respect to the values 
obtained by the Schwarzschild criterion by a quantity 
equal to 0.1 Hp, where Hp is the pressure scale height 
evaluated at the Schwarzschild boundary. 

— As in paper VIII (Meynet & Maeder 12003), 
the opacities are from Iglesias & Rogers 
H199f)|l . complemented at low temperatures 
with the molecular opacities of Alexander 
( http: / / web.physics.twsu.edu/alex/wwwdra.htmj l . 
The nuclear reaction rates are also the same as in 
paper VIII and are based on the new NACRE data 
basis (Angulo et al. ll999|l . 

— The initial composition is changed with respect to pa- 
per V: the mass fraction of hydrogen is X = 0.705, the 
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Fig. 1. Evolution of the rotational velocities for star mod- 
els of different initial masses between f 20 and 9 M© with 
account taken of anisotropic mass loss during the MS 
phase. The rapid decrease near the end of the MS evolu- 
tion is due to the bi-stability limit in the mass loss rates. 
The location of the end of the MS phase corresponds to 
the vertical decrease of the velocity. 

mass fraction of helium Y = 0.275. For the heavy ele- 
ments, we adopt the same mixture as the one used to 
compute the opacity tables for solar composition. 

The effect of rotation on the transport of the chemical 
species and the angular momentum are accounted for as 
in our papers VII and VIII. All the models were computed 
up to the end of the helium-burning phase. 

3. Evolution of the rotation 

3.1. Global view: the effects of coupling and mass loss 

The evolution of the rotation velocities at the stellar sur- 
face depends mainly on 2 factors, the internal coupling 
and the mass loss. 

1.- The coupling mechanisms transport angular mo- 
mentum in the stellar interiors. The extreme case of strong 
coupling is the classical case of solid body rotation. In this 
case when mass loss is small, the star reaches the critical 
velocity during the MS phase more or less quickly depend- 
ing on the initial rotation as shown by Sackman & Anand 
H1970I see also Langer ilQg?,! . In that case, coupling ensures 
that fJ/ilc increases when the radius of the star increases. 
Let us recall that for fi/fic to increase when the radius of 
the star increases, its suffices that ft decreases less steeply 
with the radius than flc (oc R^^^^), or expressed in an- 
other way, if f2 oc i?^" then Q/Qc increases with the ra- 
dius when a < 3/2. In the case of no couphng, i.e. of local 
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Fig. 2. Evolution of the fraction jy- of the angular veloc- 
ity to the critical angular velocity at the surface of star 
models of different initial masses between 120 and 9 Mq 
with account taken of anisotropic mass loss during the MS 
phase. 

conservation of the angular momentum, rotation becomes 
more and more subcritical (a = 2). In the present models, 
the situation is intermediate, with a moderate coupling 
due mainly to meridional circulation, which is more effi- 
cient (Meynet & Maeder |2000') than shear transport, as 
far as transport of angular momentum is concerned. 

2.- For a given degree of coupling, the mass loss rates 
play a most critical role in the evolution of the surface 
rotation. As shown by the comparison of the models at 
Z = 0.02 and Z = 0.004 (Maeder & Mevnet IM)T|l . for 
masses greater than 20 M© the models with solar compo- 
sition have velocities that decrease rather rapidly, while at 
Z = 0.004 the velocities go up. Thus, for the most massive 
stars with moderate or strong coupling, the mass loss rates 
are the main factor influencing the evolution of rotation. 

Below a mass of about 12 Mq, the mass loss rates are 
smaller and the internal coupling plays the main role in 
the evolution of the rotational velocities. This provides an 
interesting possibility of tests on the internal coupling by 
studying the differences in rotational velocities for stars 
at different distances of the ZAMS. In particular, such a 
study could allow us to test the role of magnetic coupling 
in radiative envelopes, which is now a major open question 
in stellar rotation studies f Spruit l^()()2|l . 

Figs, m and [3 show the evolution of the rotational ve- 
locities and of the fraction jy- of the angular velocity to 
the critical angular velocity at the surface of star models 
of different initial masses between 9 and 120 Mq with ac- 
count taken of anisotropic mass loss during the MS phase. 
Globally these figures are rather similar to those of pre- 
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Fig. 3. Evolution as a function of the actual mass of the 
rotation period, the rotational velocity and the ratio of 
the angular velocity to the critical value during the WR 
stage of massive stars. The continuous line corresponds 
to the so-called eWNL phase (see text), the dotted and 
short -dashed lines show the evolution during the eWNE 
and WC phases respectively. The long-dashed lines in the 
panels for the velocities show the evolution of the polar 
radius, (in this case the unity on the ordinates are R/Rq). 
Left: the WR phase of a star with an initial mass of 85 M0 
with Vini — 300 km s^^. Right: the same for an initial mass 
of 25 Mq with Vini = 300 km s~^ 



vious models (Meynet & Maeder 12000(1 . except for the 
fact that the mass loss rates by Vink et al. H2000L I200HI 
are used here. These rates meet a bi-stability limit at 
log Toff — 4.40. Below this value there is a sudden increase 
of the mass loss rates, which makes the rotation velocities 
rapidly decrease, as for the most massive stars. This is 
clearly visible for the models below 40 M0 in Figs.^andlS] 
Thus, for these models with a rapidly decreasing velocity 
of rotation, we again stress that any comparison between 
observed and predicted rotation for these domains of mass 
is really much more a test of the mass loss rates than a 
test on the internal coupling and evolution of rotation. We 
may mention that for an initial velocity Vi^i < 300 km s^^, 
the account taken of the anisotropic wind docs not play a 
major role (see Sect. 4). 



Thus, the uncertainties on the mass loss rates and bi- 
stability limits prevent us from making significant tests on 
the internal physics and coupling in the massive stars for 
models with initial masses above about 12 Mq. Mass loss 
uncertainties dominate the problem. From the rather large 
differences between the values of the mass loss rates pub- 
lished over recent years, we know that these uncertainties 
are still great. Fortunately, below 12 Mq the comparisons 
may be meaningful and may reveal whether our models 
have the right kind of internal coupling. 

3.2. Evolution of the rotational velocities during the 
WR stages 

Fig. 131 (left) shows the evolution of the rotation period 
P, of the rotation velocities v and of the fraction ^ of 
the angular velocity to the critical angular velocity at the 
surface during the WR stages of a star model with an 
initial mass of 85 Mq with Vini = 300 km s^^ typical of a 
star which has both a WN and WC phase. During the WN 
stage of the 85 Mq model, until the actual mass of the star 
reaches a value near 23 Mq , the angular velocity decreases 
by a factor of ~ 5, correspondingly the period P = 
grows by the same factor and the rotation velocity v = 
flRo decreases much more since the equatorial radius Re 

also decreases during this phase (by a factor of 2). The 
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also decreases very much, because 



fraction M- 

(GM)5 

the sensitivity to the decrease of the radius is larger than 
that to the mass (which decreases by a factor of ^ 2.7). 

We adopt here the definition of the WN phases given 
by Foellmi et al. I|2003|l . putting a prefix "e" when the class 
is assigned by composition properties of the evolutionary 
models. Thus, eWNE means WN stars without hydrogen, 
while eWNL means WN stars with hydrogen, (see also 
Sect 5 below). When the H-content goes to zero (i.e. just 
before the star enters the eWNE stage) , we see a dramatic 
increase of (or decrease of P) , due to the strong contrac- 
tion of the surface layer because the opacity of the surface 
layers becomes much smaller. This transition is fast, thus 
the transfer of angular momentum by meridional circu- 
lation is limited and the evolution of the rotation at the 
surface is dominated by the local conservation of angular 
momentum, which explains the rapid and large increase of 
ri. However, we see that despite the strong decrease of P, 
the velocity v and the ratio ^ do not change very much, 
since the radius is much smaller in the eWNE stage. When 
the star becomes a WC star with an actual mass of ~ 17.3 
M0, there is very little change in the velocity, since rota- 
tion is already very low and the radius no longer decreases 
significantly. 

Fig. El (right) shows the same for a star of 25 Mq with 
Wini — 300 km s^^, which does not enter the eWNE and 
WC stages. We see that P decreases (or il increases) dur- 
ing the WN stage. This is because the star model moves 
from log Tcff 4.29 to 4.70. The decrease of the radius makes 
the star rotate faster, slightly overwhelming the mass loss 
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Fig. 4. Same as Fig.Elfor 60 Mq stellar models. Left: the 
WR phase of a star with an initial mass of 60 M© and 
Willi — 300 km s^^. Right: for an initial mass of 60 Mq 
with iiini = 500 km s~^. In the e-version, the red color 
corresponds to the phase where H is present (eWNL), the 
black to the phase of an He-star (eWNE) and the blue to 
the WC phase. 

effects from 13.65 to 11.33 Mq. The increase of fl and the 
decrease of R nearly compensate for each other so that v 
does not vary much during most of the WN phase, except 
at the very end. There the increase of Q is more rapid 
than the decrease of the radius and an increase in the ve- 
locity is obtained. The reason for this behaviour is that 
in this model the stellar winds uncovers the layers, where 
the ri-gradient is steep, and we start seeing layers which 
were deep inside and rotating rapidly. 

Globally, we note that the rotational velocities of WR 
stars are modest in the 2 cases examined above, despite 
the very different masses and evolution stages. This is 
firstly the result of the strong mass loss necessary to pro- 
duce the WR stars and also of the mass loss during the 
WR stages. We note however that, for a given initial veloc- 
ity, the velocities obtained at the end of the He-burning 
phase for the WR stars decrease when the initial mass 
increases (see Table This results from the simple fact 
that the higher the initial mass, the larger the quantity 
of mass (and thus of angular momentum) removed by the 
stellar winds. 



Figs. 01 shows the evolution of rotation properties for 
a 60 Mq model with Vini = 300 km (left) and with 
Wini = 500 km s~^ (right). On the whole, these two cases 
are not very different from the case of the 85 Mq model 
shown above. In Fig.0](left), we firstly see an increase of P 
and a corresponding decrease of and of ^ ; this decrease 
of the rotation is due, as before, to the high mass loss rates. 
Then, when the actual mass is about 38 Mq, 57 increases 
again due a decrease of the stellar radius as the star is 
moving to the left in the HR diagram (Fig.CJ. The velocity 
V does not change very much since the stellar radius is 
decreasing. When the actual mass of the star approaches 
30 M0, at an age of 4.307- 10^ yr, = 0.10 and the outer 
stellar layers become more transparent, so that the star 
brightens rapidly above logL/L© = 6.0 for a short time 
until an age of 4.33 • 10^ yr. The corresponding variations 
in radius produce a sharp oscillation in Q and velocity. 
Then at an age of 4.38 • 10^ yr. until 4.44 • 10^ yr., the star 
becomes an He-star (eWNE stage). The angular velocity 
remains nearly constant (cf. the evolution of the period) 
while the velocity and Q/flc present some more variations 
linked to changes in the radius. Then, the star enters the 
WC stage, the mass loss rates are so large that despite 
the decrease of the radius, the rotation remains small. As 
was the case before, the presence of steep gradients of J7 
at the surface, may result in some cases in an increase of 
fl when the star is peeled off by the stellar winds. 

The case of 60 Mq with Vi^i — 500 km s~^ (right) con- 
tains features similar to the previous cases. The beginning 
of the curve, with the first bump of P, behaves like the 
case of Wini — 300 km of 60 M©. After the dip of P at 
an age of 4.08-10^ yr (actual mass of ~ 37 M©), the further 
evolution of P is like the case of an 85 M© star, marked 
by a large decrease of the angular velocity Q. Later, when 
the star becomes an He-star at an age 4.50 • 10^ yr (ac- 
tual mass 24 M©), we see a sharp increase in il, in the 
velocity and in the ratio il/ilc- During the WC phase, the 
velocity remains between 30 and 40 km s^^, surprisingly, 
this is a factor ~ 2 below the velocity reached during this 
phase by the 300 km s^^ stellar model. The lower value 
obtained for the initially faster rotating star results from 
the greater quantity of mass and thus of angular momen- 
tum lost during the previous evolutionary stages. 

The main conclusion of this Section is that the pre- 
dicted surface rotation velocities of WR stars are relatively 
low, of the order of 50 km s^^. Also, in the stellar interiors, 
rotation is modest with ratios of the angular velocity 
at the center to that at the surface which are only of the 
order of a few unities. This is the result of the very strong 
mass loss peeled off the star to make it a WR star. It is 
to be examined whether the further evolution and strong 
central contraction to the pre-supernova stage is able to 
accelerate enough to account for the observed rotation pe- 
riods of pulsars and maybe also, for rotations fast enough, 
for enabling the coUapsar model to work fWooslev l2003|l . 
A word of caution is necessary, since there are many paths 
leading to WR stars and we cannot exclude that some, 
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Fig. 5. Evolutionary tracks for non-rotating (dotted lines) and rotating (continuous lines) models for a metallicity Z 
= 0.020. The rotating models have an initial velocity Wini of 300 km s^^, which corresponds to an average velocity 
during the MS phase of about 180 to 240 km s"^ (see TableCJ. 



possibly rare paths, may lead to WR stars with very high 
rotation. 

4. Evolutionary tracks, masses and lifetimes 

Fig. |S1 presents the tracks for non-rotating and rotating 
evolutionary models for the whole mass range considered. 
As a result of the change of the opacities (slightly en- 
hanced) and of the initial abundances {Y = 0.275, instead 
of 0.30), the non-rotating tracks are shifted to lower ef- 
fective temperatures with respect to paper V. The tracks 
also extend towards lower effective temperatures during 
the MS phase. This is a well known consequence of over- 
shooting. Overshooting and core extension by rotation are 
also generally responsible for shorter extensions of the blue 
loops in the He-burning phase. This is what we observe in 
Fig. 13 the extension of the blue loops is further reduced 
by rotation in the case of the 9 Mq model, and the blue 
loops are even suppressed in the rotating model of 12 M©. 

The present non-rotating 12 M©, with overshooting, 
presents a blue loop, while the non-rotating model of pa- 



per V, without overshooting, shows no blue loop. This is at 
first sight surprising because, as said above, more massive 
He-cores tend to reduce the extension of the blue loops, 
if not suppressing them (Lauterborn et al. 119711 see also 
the discussion in Maeder & Me vnet [1989(1 . As explained in 
paper V, this particular situation, at the limit where the 
loops generally disappear, results from the higher luminos- 
ity of the model having the more massive He-cores (either 
as a result of overshooting or as a result of rotation) . The 
higher luminosity implies that the outer envelope is more 
extended, and is thus characterized by lower temperatures 
and higher opacities at a given mass coordinate. As a con- 
sequence, in the model with higher luminosity, during the 
first dredge-up, the outer convective zone proceeds much 
more deeply in mass, preventing the He-core to grow too 
much during the further evolutionary phases and enabling 
the apparition of a blue loop. 

The effects of rotation on the tracks have already 
been discussed in previous works fHeeer [r998l Meynet & 
Maeder r2()()()|l . thus we shall be very brief here. Let us just 
list the main effects: 
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Table 1. Properties of the stellar models at the end of the H-burning phase and at the end of the He-burning phase. 
The masses are in solar mass, the velocities in km s^^, the lifetimes in million years and the abundances in mass 
fraction. 
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^ The symbol I after the velocity indicates that the model was computed assuming isotropic stellar winds. 



^ This model was computed until the mass fraction of hydrogen at the center is 0.000348. 

^ The minimum initial mass star of O— type star is 15.9 Mq for non— rotating models and 17 Mq for the rotating stellar models. 



— On and near the ZAMS, the atmospheric effects of ro- 
tation produce a small shift of the tracks towards lower 
Toff. The effective temperature considered here is an 
appropriate averaged Tgff over the stellar surface (see 
Meynet & Maeder [TM7|l . 

— Tracks with rotation become more luminous than for 
no rotation. 

— For stars with initial masses inferior to ^ 30 M0, the 
MS bandwidth is enlarged by rotation. This is a con- 
sequence of the larger He-core produced during the 
H-burning phase in rotating models. As a numerical 
example, in the 20 Mq model with Vjni = 300 km s~^, 
the He-core at the end of the MS is about 38% more 
massive than in its non-rotating counterpart. In the 
models of paper V the enhancement due to rotation 
was only 20%. Likely most of the difference results 
from the inclusion here of a small overshooting. In the 
present models, the core enhancement due to rotation 
adds its effects to an already more extended core, hav- 



ing a higher central temperature, and thus this larger 
core reacts more strongly to further supplies of fuel by 
rotational diffusion. 
— For initial masses superior to about 50 M0 , the rotat- 
ing tracks show drastically reduced MS bandwidths. 
This was already the case in Meynet & Maeder H2000fl . 
For instance, the reddest point reached at the end of 
the MS phase by the 60 M0 non-rotating stellar model 
is log Tcff ^ 4.15 , while for the Vini — 300 km s^^ 
model, it is ~ 4.6 , due to the fact that the surface He- 
enhancement leads to a bluer track, as a consequence 
of the reduced opacity. In this mass range, the rotat- 
ing stellar models enter into the WR phase during the 
H-burning phase (see Sect. 5). The MS phase, iden- 
tified here with the OV-type star phase, is thus con- 
siderably reduced. Since the stars, here as in any mass 
range, have different initial rotational velocities, one 
expects to find in the HR diagram stars correspond- 
ing to a variety of tracks with various initial velocities. 
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Fig. 6. Relations between the final mass versus the ini- 
tial mass for solar metallicity models. The cases with and 
without rotation are indicated. The line with slope one, la- 
beled by M = 0, would correspond to a case without mass 
loss. The relation obtained from the models of Schaller et 
al. I|1992|l is also shown. 
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Particularly, in the upper hot part of the HR diagram, 
for log L/Lq superior to about 5.8-5.9 and log Teg 
between 4.7 and 4.2, one may find in the same area 
of the HR diagram stars in very different evolutionary 
stages, which may also correspond to different types 
of massive stars: MS 0-type stars. Of stars, transition 
stars, blue supergiants and WR stars. 

For higher rotation velocities, the actual masses at the 
end of both the H- and Hc-burning phases become smaller 
(cf. Tables ^1. For an initial velocity of 300 km s^^ and 
for stars with masses below ~ 60 M0, the actual masses 
at the end of the core H-burning phase are decreased by 
less than 11% with respect to the non-rotating models. 
For higher initial mass models, the decreases are much 
more important, because these stars enter the WR phase 
already at an early stage during the H-burning phase (see 
Sect. 5). 

The evolutionary masses of a star can be determined 
by searching the mass of the evolutionary track passing 
through its observed position in the HR diagram. As can 
be seen from Fig. [3 the mass determined from the non- 
rotating evolutionary track is always greater than the one 
deduced from rotating models. Thus the use of non rotat- 
ing tracks tends to overestimate the mass, and this might 
be a cause of a large part of the long-standing mass dis- 
crepancy problem (see e.g. Herrero et al. 12000(1 . 

Fig.|n|shows the final masses obtained for different ini- 
tial masses and velocities. One notes that the new mass 



Fig. 7. The two upper panels show the evolutionary tracks 
of a non-rotating and a rotating 60 M© stellar model. 
The bold part of the track corresponds to the 0-type star 
Main-Sequence phase; the dotted part shows the inter- 
mediate phase, when present, between the 0-type star 
phase and the WR phase; the thin continuous line is the 
track during the WR phase. The two lower panels show 
the evolution as a function of time of the total mass and 
of the masses of the convective cores during the H- and 
the He-burning phases. 

loss rates used in the present computations are smaller 
than those used in the stellar grids of Schaller et al. H1992fl . 
For stars with initial masses above 40 M0 , rotating models 
produce final masses between 11.3 and 14.7. For compar- 
ison, non-rotating models end with final masses between 
14.1-17.3 Mq. Interestingly, the average mass of WC stars, 
estimated on the basis of the 6 WC stars belonging to 
double-line spectroscopic binaries, is 12 ± 3 M© (van der 
Hucht .2001f) . nearer to the rotating model results than the 
non-rotating ones. However we cannot discard the possi- 
bility that for these stars mass transfer in a Roche lobe 
overflow may have occurred, in which case this observation 
cannot be used to constrain single star models. The fact 
that the models with and without rotation near 50 Mq 
lead to about the same final masses, despite the increased 
mass loss rates for rotating models, is due to the differ- 
ences in the evolutionary stages followed by each model 
(see Sect. 5). 
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Fig. 8. Evolutionary tracks for a non-rotating (light 
dotted-line) and a rotating 120 Mq stellar model at solar 
metallicity. The circles along the tracks indicate the end 
of the eWNL phase, the squares the end of the eWNE 
phase, and the triangles, the beginning of the WC phase. 
The WN/WC stars lie along the portion of the tracks be- 
tween the squares and the triangles. The crosses indicate 
the end of the central H-burning phase. 



is much greater and thus rotation predicts larger number 
ratios of WR to 0-type stars than non-rotating stellar 
models. 




Fig. 9. Lifetimes of WR stars of different initial mass for 
non-rotating stellar models. The durations of the WR sub- 
phases are also indicated. There is no transition WN/WC 
phase. 



Table ^ presents some properties of the models. 
Columns 1 and 2 give the initial mass and the initial veloc- 
ity Willi respectively. The mean equatorial velocity v during 
the MS phase is indicated in column 3. This quantity is 
defined as in paper V. The H-burning lifetimes tn, the 
lifetimes to as an 0-type star on the MS (we assumed 
that 0-type stars have an effective temperature superior 
to about 31 500 K as suggested by the new effective tem- 
perature scale given by Martins et al. I2002|l . the masses 
M, the equatorial velocities v, the helium surface abun- 
dance Ys and the surface ratios (in mass) N/C and N/0 
at the end of the H-burning phase are given in columns 4 
to 10. The columns 11 to 16 present some characteristics 
of the stellar models at the end of the He-burning phase, 
tHe is the He-burning lifetime. The letter I attached to 
the value of the initial velocity (column 2) indicates that 
the model was calculated assuming isotropic stellar winds. 
These models were calculated until the end of the MS 
phase, except in the case of the 60 M© model for which 
the evolution was pursued until the end of the core He- 
burning phase. 

From Tablcnone sees that for Z = 0.020 the lifetimes 
are increased by about 15-25% when the initial rotational 
velocity on the MS increases from to 300 km s^^. Similar 
increases were found in paper V. The lifetimes of 0-type 
stars are increased by similar amounts. As we shall see 
below, the increase due to rotation of the WR lifetimes 



Except for the models with initial masses superior to 
85 Mq, the He-burning lifetimes are decreased by rota- 
tion. Thus the ratio of the He-burning lifetimes to the 
H-burning lifetimes are in general decreased by rotation. 
Typically for the 20 Mq model, this ratio passes from 
0.11 to 0.08 when the initial velocity considered increases 
from to 300 km s~^. The particular behaviour of the ex- 
tremely massive stars results from the lower luminosities 
in the WC stage, which favour longer He-lifetimes. 

The anisotropy of the wind induced by rotation has 
little impact when the initial velocity is equal or inferior 
to 300 km s~^. This can be seen by comparing the data 
presented in Tabled for the two 60 Mq models with Vini 
= 300 km computed with and without the account of 
the anisotropy. They are indeed very similar. On the other 
hand, as has been shown by Maeder H2()()2|l . for higher 
initial velocities, in the mass range between 20 and 60 
Mq, the anisotropics of the mass loss may significantly 
influence the evolution of the stellar velocities and may 
lead many stars to break-up. 

5. Rotation and the properties of Wolf-Rayet 
stars 

Here we consider the following question: what are the ef- 
fects of rotation on the evolution of massive single stars 
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Fig. 10. Lifetimes for rotating stellar models in the WR 
subphases. 

into the Wolf-Rayet phase ? This subject has been dis- 
cussed by Maeder fT987), Fliegner and Langcr (1995), 
Maeder & Meynet f2000b) and Meynet ( 2000 ). We present 
here for the first time an extensive grid of models, account- 
ing for the effects of the wind anisotropy, going through 
the WR phase. It is important to cover the whole mass 
range of star progenitors of WR stars in order to quanti- 
tatively assess the importance of rotation in the interpre- 
tation of the WR number statistics (see Sect. 6 below). 

Let us reconsider the criteria which have been chosen 
to decide when a stellar model enters into the WR phase. 
Ideally, of course, the physics of the models should decide 
when the star is a WR star. However our poor knowledge 
of the atmospheric physics involved, as well as the com- 
plexity of models coupling stellar interiors to the winds, 
are such that this approach is not yet possible. Instead, it 
is necessary to adopt some empirical criteria for deciding 
when a star enters the WR phase. In this work the criteria 
are the following: the star is considered to be a WR star 
when logTcff > 4.0 and the mass fraction of hydrogen at 
the surface is Xs < 0.4. Reasonable changes to these val- 
ues (for instance adopting Xs < 0.3 instead of 0.4) do not 
affect the results significantly. 

The correspondence between the WR-subclasses and 
the surface abundances is established as follows: the stars 
enter the WR phase as WN stars i.e. as stars presenting 
strong He and N enrichments at their surface and strong 
depletion in carbon and oxygen as expected from CNO 
processed material. We call eWNL stars all the WN stars 
showing hydrogen at their surface ("e" stands for evo- 
lutionary, as recommended by Foellmi et al. 12003(1 . The 
eWNE stars, on the other hand, are WN stars showing 
no hydrogen at their surface. The transitional WN/WC 



phase is characterized by the simultaneous presence at the 
surface of both H- and He-burning products (typically N 
and C or Ne). Here we suppose that this transition phase 
begins when the mass fraction of carbon at the surface 
becomes superior to 10% of the mass fraction of nitro- 
gen. The phase is supposed to end when the mass frac- 
tion of nitrogen at the surface becomes inferior to 10% of 
the mass fraction of carbon. The WC phase corresponds 
to the stage where He-burning products appear at the 
surface. This phase begins when the transition WN/WC 
phase ends. 



Wolf-Rayet lifetimes 
of 60 Mq models 




200 



400 



[km s-'] 



Fig. 11. Durations of the WR phases and of the WR sub- 
phases for stars of 60 M© with various initial velocities. 

Figure [71 shows the evolutionary tracks and the evo- 
lution of the structure in a non-rotating and a rotating 
60 Mq model. The rotating model has a time-averaged 
equatorial velocity on the Main Sequence (MS) of about 
190 km s^^, not far from the mean equatorial velocity ob- 
served for 0-type stars. Inspection of Fig. [7\ shows the 
following features: 

— The non-rotating 60 M© star goes through an inter- 
mediate Luminous Blue Variable (LBV) phase before 
becoming a WR star, while the rotating model skips 
the LBV phase and goes directly from the 0-type star 
phase into the WR phase. One also notes that dur- 
ing the 0-type star phase, the rotating track is bluer 
than its non-rotating counterpart. This is due to the 
diffusion of helium into the outer radiative zone (see 
Heger & Langer 2000; Meynet & Maeder EOT!) . This 
fact may slightly change the number of 0-type stars 
deduced from the number of ionizing photons, a tech- 
nique used to estimate the number of 0-type stars 
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from the spectral analysis of the integrated light of 
distant starburst regions (e.g. SchaererfTWBI. 

— There is no luminosity difference during the WC phase 
between the rotating and the non-rotating models. 
Since these stars follow a mass-luminosity relation 
(Maeder [TMTl Schaerer & Maeder HM^ . this im- 
plies that the masses of these two stars at this stage 
are about equal (see Fig. At first sight this may 
be surprising. One would expect the rotating star 
to lose more mass than its non-rotating counterpart. 
However, in this particular case there are some com- 
pensating effects. Indeed, the non-rotating model en- 
ters into the WR phase later, but becomes redder 
during the MS phase and then becomes a LBV star. 
During these two last phases, the star undergoes strong 
mass loss (see the left panels of Fig. [T)). The rotating 
model, on the other hand, enters the WR phase at an 
earlier stage, reaching at the end an identical mass to 
the non-rotating model. In general, however, the final 
masses obtained from the rotating models are smaller 
than those obtained from the non-rotating ones (see 
Fig. O, leading to lower luminosities during the WC 
phase (see Fig. ISJ. 

— From the lower panels of Fig. |7| we see that the non- 
rotating star enters the WR phase at the beginning 
of the core He-burning phase with a mass of about 
27 Mq. Nearly the whole H-rich envelope has been 
removed by stellar winds during the previous phases 
(more precisely at the end of the MS phase and during 
the LBV phase). In this case, the main mechanism for 
making a WR star is mass loss by stellar winds. For the 
rotating model, on the other hand, the entry into the 
WR phase occurs at an earlier stage (although not at a 
younger age!), while the star is still burning hydrogen 
in its core. The mass of the star at this stage is about 
45 Mq and a large H-rich envelope is still surrounding 
the convective core. The main effect making a WR 
star in this second case is rotationally induced mixing 
("Maeder [TWI Fliegner & Langer llMSI Meynet 110001) • 

In Fig. IHl the evolutionary tracks for a non-rotating 
and a rotating 120 Mq stellar model are shown. The be- 
ginning and the end of the various WR subphases are indi- 
cated along the tracks. The range of luminosities spanned 
by the rotating eWNL stars is considerably increased. This 
directly results from the point underlined just above, i.e. 
a rotating star enters the WR phase still with a large H- 
rich envelope. This considerably increases the duration of 
the eWNL phase and the star will cover during this stage 
a greater range of masses and luminosities. The rotating 
models also enter the further WR subclasses with smaller 
actual masses and lower luminosities. From Fig.|Hl one sees 
that the rotating eWNE, WN/WC and WC stars have lu- 
minosities decreased by about a factor of two. 

The differences between the behaviours of the rotating 
and non-rotating models have important consequences for 
the duration of the WR phase as a whole and for the 



lifetimes spent in the different WR subtypes, as can be 
appreciated from Figs. 151 and 1101 

— Firstly, one notes that the WR lifetimes are enhanced. 
Typically for the 60 Mq model, the enhancement by 
rotation amounts to nearly a factor of 2. 

— Second, the eWNL phase is considerably lengthened 
as explained just above. The duration of this phase in- 
creases with the initial velocity as can be seen from 
Fig. When the initial velocity passes from 300 to 
500 km s^^ (or when the time-averaged equatorial ve- 
locity during the 0-type phase goes from 190 to 330 
km s~^), the eWNL phase duration is increased by a 
factor of about 1.7. The durations of the eWNE and 
WC phases are, on the other hand, much less affected. 

— Thirdly, in the rotating stellar model a new phase of 
modest, but of non-negligible duration, appears: the 
so-called transitional WN/WC phase. This phase is 
characterized by the simultaneous presence at the sur- 
face of both H- and He-burning products. The reason 
for this is the shallower chemical gradients that build 
up inside the rotating models. These shallower gradi- 
ents inside the stars also produce a smoother evolu- 
tion of the surface abundances as a function of time 
(or as a function of the remaining mass, see Figs. ITHl 
and in^ . For a transitional WN/WC phase to occur, it 
is necessary to have — for a sufficiently long period — 
both a He-burning core and a CNO-enriched enve- 
lope. In the highest mass stars, mass loss removes too 
rapidly the CNO-enriched envelope to allow a tran- 
sitional WN/WC phase to appear. In the low mass 
range, the time spent in the WR phase is too short and 
the H-rich envelope too extended to allow He-burning 
products to diffuse up to the surface. Consequently, the 
transitional WN/WC phase only appears in the mass 
range between ^ 30 and 60 M0 for Wini — 300 km s~^. 

— Finally, the minimum mass for a star to become a WR 
star (through the single star channel) is lowered by 
rotation. In the present case, the minimum mass is 
reduced from about 37 M© for Vjni = km s~^ to 
about 22 Mq for Vini = 300 km s^^. 

Of course, this applies to solar metallicity. At very dif- 
ferent Z, the effects of rotation on the various subphases 
may be different. It is interesting to compare the present 
WR-lifetimes with those we obtained with our previous 
non-rotating stellar models (Meynet et al. 119941 Maeder 
and Meynet Unnil). This is done in Fig. [El We see that 
the present non-rotating models follow more or less the 
behaviour we obtained with our normal mass loss rate 
models of 1994, although, due to smaller mass loss rates in 
the present grid, the WR lifetimes are in general reduced 
and the minimum mass for WR star formation from sin- 
gle stars is enhanced. Our rotating models, on the other 
hand, reproduce qualitatively the changes brought by an 
enhancement of the mass loss rates: an increase of the WR 
lifetimes and a decrease of the minimum mass for a single 
star to go through a WR phase. 
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From the above considerations, one can conclude that, 
compared to non-rotating models, stellar models includ- 
ing rotation will predict higher values of: (a) the numbers 
of WR stars relative to 0-type stars, (b) the relative num- 
ber of WN to WC stars, and (c) the relative number of 
transition WN/WC stars to WR stars. As shown in the 
next Section, the predictions of the rotating models are 
in better agreement with the observed number ratios at 
solar metallicity. 
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Table 2. Lifetimes of Wolf-Rayet star from various initial 
masses and velocities and predictions for different number 
ratios (see text). The velocities are in km s^^ and the 
lifetimes in Myr. 
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Fig. 12. Lifetimes of Wolf-Rayet stars from various ini- 
tial masses at solar metallicity. The continuous line shows 
the results from the present rotating models. The WR 
lifetimes of the non-rotating models from the present 
work are shown by the dotted line. The short- and long- 
dashed curves show the results obtained from the models 
of Schaller et al. ijlM^ and Meynet et al. ifDMjl with 
normal (1 X Al) and enhanced mass loss rates (2 X M) 
respectively. 



6. Comparison with the observed properties of 
Wolf-Rayet star populations 

Following Maeder (fTMT]l and Maeder & Meynet lfBM| . 
one can easily estimate the theoretical number ratio of 
WR to 0-type stars in a region of constant star forma- 
tion. This ratio is simply given by the ratios of the mean 
lifetime of a WR star to that of an OV-type star. The av- 
erages are computed using an initial mass function as the 
weighting factor over the respective mass intervals consid- 
ered. Assuming a Salpeter Initial Mass Function slope of 
1.35, considering the 0-type and WR star lifetimes given 
in tables ^ and 12 we obtain the predicted ratios given at 
the bottom of this last table. 



Comparisons with observed number ratios are shown 
in Figs. ll3l andll4l Obviously, unless most of the WR stars 
at solar metallicity are formed through Roche lobe over- 
flow in a close binary system, which seems quite unrealistic 
(see the discussion in Maeder I199(jl see also Foellmi et al. 
I2U(J3|I , models with rotation are in much better agreement 
with observations. Of course a more refined estimate of 
the effects of rotation will require account of the initial 
distribution of the rotational velocities. However, since 
the initial velocity chosen here corresponds to a mean ro- 
tational velocity on the Main Sequence almost equal to 
the observed mean, we can hope that the present result 
is not too far from what one would obtain with a more 
complete analysis. The conclusion is that at solar metal- 
licity, rotating models are in much better agreement with 
the observed statistics of WR stars than the non-rotating 
models. Interestingly a similar conclusion could be reached 
based on the results of our paper V, and this shows that 
the present results depend weakly on overshooting or on 
the inclusion of the effects of horizontal turbulence in the 
coefficient of shear diffusion. 

In Fig. nil the predicted values for the WC/WN num- 
ber ratios are compared with observations. One sees that 
the WC/WN ratio increases with the metallicity along 
a relatively well defined relation. The observed point for 
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Fig. 13. Variation of the number ratios of Wolf-Rayet 
stars to 0-type stars as a function of the metaUicity. The 
observed points are taken as in Maeder & Meynet (jlPMj) . 
The dotted and continuous lines show the predictions of 
the models of Schaller et al. 1)1992(1 and Meynet et al. 
H1994|l with normal and enhanced mass loss rates respec- 
tively. The black square and triangle show the predictions 
of the present non-rotating and rotating stellar models 
respectively. 

the solar neighborhood is however well above the general 
trend. According to Massey (2003) this may result from 
an underestimate of the number of WN stars. 

Comparing with the theoretical predictions at solar 
metallicity, one obtains that the non-rotating models pre- 
dict too few WN stars, while the rotating models give a 
value in the observed trend. This confirms the conclusion 
already reached above on the basis of the WR/0 num- 
ber ratio: rotating models much better account for the 
characteristics of the WR populations than non-rotating 
models. 

There are 10 WN/WC stars listed in the Vllth cat- 
alogue of galactic WR stars of van der Hucht (2001). 
Compared to the total number of known galactic WR 
stars (226) this represents a fraction of 4.4%. The present 
models predict for this fraction a value equal to 4% in 
good agreement with the observed ratio. We recall that 
models without mild mixing processes in radiative zones 
are unable to predict correctly the number of transition 
WN/WC stars. This was first shown by Langer H1991|l . 

In future work, it has to be checked by further compu- 
tations at different metallicities that the rotating models, 
with reduced mass loss rates due to clumping, will account 
for the observed variation of the number ratios WR/0 and 
WC/WN with metaUicity. 

7. Surface abundances in WR stars 

We have already presented and discussed the effects of ro- 
tation on the surface abundance of massive stars at solar 
composition (Meynet & Maeder r2()()()|l . The main point is 



the enhancement of He and N during the MS phase due to 
the internal transport by shear mixing. However, we note 
that the enhancements in He and CNO products found 
here are in general smaller that those found in Meynet 
& Maeder The situation is illustrated in Fig. [13 

which compares the evolution of the N/C ratios for mod- 
els with masses between 9 and 40 Mq of the above quoted 
paper and of the present work. We see that a difference in 
the enhancement of N/C is created during MS evolution 
and that the present ratios are smaller on average by a 
factor of about 1.5 for masses below 25 M0. This is due 
to the fact that in the present models we are accounting 
for the effects of the horizontal turbulence on the shear 
mixing according to the study by Talon & Zahn H1997|) . 
The horizontal diffusion substantially reduces the effects 
of the shear mixing. In the higher mass range, the dif- 
ference between the two sets of models is smaller, as can 
be seen by comparing the two 40 M© tracks in Fig. 
In this domain of masses, the large effects of mass loss 
by stellar winds are added to those of rotational mixing. 
Close comparisons with observations will be undertaken 
on the basis of these two sets of models to see whether the 
treatment by Talon & Zahn leads to a better agreement or 
not. However, the problem is intricate both observation- 
ally, because the scatter of the observations is large, and 
also theoretically because other effects may enter into the 
game, such as magnetic fields (Spruit ' 2002|) or a higher 
horizontal turbulence. Let us also recall that the enhance- 
ments in He and N are much larger at lower metallicities 
Z (Maeder & Mevnet [2()()l|l . because stars with lower Z 
have in general steeper fi-gradients, which favour internal 
mixing. 

Fig. 1 161 shows the evolution of the surface abundances 
in a rotating and a non-rotating model with initial masses 
of 60 Mq . The progressive changes of the abundances of 
CNO elements between the initial cosmic values and the 
values of the nuclear equilibrium of the CNO cycle are 
much smoother for the rotating model than for the case 
without rotation. This is due to the internal mixing which 
makes flatter internal chemical gradients for rotating mod- 
els. The rotational mixing also makes the change of abun- 
dances to occur much earlier in the peeling-off process. 
This is also true for the changes of H and He. Another 
significant difference for the 60 Mq model with rotation is 
the absence of an LBV phase as shown in Sect. 5. Without 
rotation, most of the change of CNO abundances occurs 
in the LBV phase, while in the case with rotation, the 
transition from cosmic to CNO equilibrium abundances 
occurs partly during the MS phase and partly during the 
eWNL phase which immediately follows it. We note that 
the nuclear equilibrium CNO values are essentially model 
independent as already stressed a long time ago (Smith & 
Maeder Il991|l . This is true whether H is still present or 
not. 

We see that at the end of the assumed eWNE phase, 
the transition to the WC abundances is very sharp in 
the absence of rotational mixing. This is because there 
is a strong chemical discontinuity in the classical models 
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Fig. 15. Evolution of the N/C ratios (in number) as a 
function of the effective temperature at the surface of vari- 
ous rotating stellar models. The N/C ratios are normalized 
to their initial values. The continuous lines show the re- 
sults for the present rotating models (uini — 300 km s~^), 
the dashed-lines show the evolutionary tracks for some 
models of paper V. For purpose of clarity, only a part of 
the 40 Mq tracks are shown. 



due to the fact that the convective core is usually grow- 
ing during some part of the He-burning phase (Maeder 
1983). In the rotating case, the transitions are smoother, 
so that there are some stars observed with abundances 
in the transition regime, with simultaneously some ^^C 
and some ^''N present. As stated above, these correspond 
to the so-called WN/WC stars (Conti & Massev IT^ 
Crowther et al. 11995(1 . These transition stars are likely to 
also have some ^^Ne excess. Since the attribution of spec- 
tral types is a complex matter, it may even be that some 
of the stars in the transition stage are given a spectral 
type eWNE or WC, thus we might well have a situation 
where a WN star would have some ^^Ne excess or a WC 
star would still have some ^"^N present. 

Fig. El compares the predicted relation between the 
C/N and C/He ratios with some recent observations by 
Crowther et al. |J995). During the eWNL phase, C/He 
and C/N decreases, but remain around the CNO equi- 
librium value equal to 0.015. Then during the transi- 
tional WN/WC phase, there is a huge increase of car- 
bon, while nitrogen and helium remain constant. This is 
the reason why the lower part of the curve shows that 
when C/N changes by a factor 10'^, C/He changes by the 
same amount. Then, when the star evolves towards the 
WC stage, C/He changes only slowly, while C/N grows 
very rapidly, since N quickly disappears. The comparison 
with the observations by Crowther et al. 1)1995(1 shows a 
good agreement. However this agreement may be some- 
what misleading here, in the sense that a 60 Mq non- 
rotating model would follow a similar track in the C/N 
versus C/He plane. The good correspondence here is not 
a real constraint on the models. What is more constrain- 
ing is the expected number of stars in the transitional 
WN/WC phase. As seen in Sect. 5, the rotating models in 
this respect are much better than the non-rotating ones. 

In the WC stage, the abundances of chemical elements 
are very much model dependent as was first shown by 
Smith & Maeder 1(1991)) . This is because in the WC stars, 
we see the products of partial He-burning. At the entry 
in the WC phase, the ^^C and ^^O abundances are lower, 
and that of He is higher in the models with rotation, due 
to their smooth chemical gradients. Also, the fraction of 
the WC phase spent with lower C/He and 0/He ratios 
is longer in models with rotation. However, rotation does 
not affect the high ^^Ne abundance during the WC phase. 
This is a consequence of the fact that most of this ^^Ne 
results from the transformation of the ^''N produced by 
the CNO cycle in the previous H-burning core. The value 
of the ^''N and therefore that of the ■^^Ne is fixed by the 
characteristics of the CNO at equilibrium, which in turn 
depends on the nuclear physics and not on the peculari- 
ties of the stellar models. It is interesting to mention here 
that the high overabundance of ^^Ne at the surface of the 
WC star predicted by the models is confirmed by the ob- 
servations (Willis 1999, Dessart et al. 12000)1 . In the core, 
the abundance of ^^Ne will stay at this equilibrium value 
until eventually it is turned into 25,26]yjg ^[^^^ production 
of neutrons (through ^^Ne(Q!,n)^^Mg) and of s-elements. 
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Fig. 16. Evolution as a function of the actual mass of the 
abundances (in mass fraction) at the surface of a non- 
rotating (upper panel) and a rotating (lower panel) 60 
M0 stellar model. 

Fig. 1181 shows the same for models of 40 Mq with and 
without rotation. The enhancements of He and N only ap- 
pear in the blue and yellow supergiant stage in the case 
without rotation, while they occur during the MS and blue 
supergiant phase when rotation is included; also in this 
case the changes are more progressive. The non-rotating 
model has almost no WNL phase and does not enter the 
WC stage. The model with rotation is very different, it 
has an extended WNL phase and then a remarkably large 
transition phase, where ^^N is present with a significant 
abundance, while ^^C, ^^O and ^^Ne are also there. Most 
of this transition phase is likely spent in the transition 
WN/WC phase, but probably not the whole, so that some 
peculiar WN stars with products of He-burning are pos- 
sible. Also some WC stars with an unusual occurrence of 
^*N are very likely. Globally, we see that the transition 
phase where both characteristics of H- and He-burning 
phases are present is longer for stars of lower masses (see 
Sect. 5). 

Fig. 1191 shows an interesting diagram for the study of 
so-called "slash stars" of type Ofpe/WN, of the LBV and 
of the WN stars. It presents the evolution of the surface 
H-content as a function of the luminosity. The tracks go 
downwards in this diagram; first there is an initial bright- 
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Fig. 18. Same as Fig.^jfor 40 Mq stellar models. 



ening due to MS evolution, then the luminosity is about 
constant, except for < 0.20 if the mass loss is very high. 
Thus for stars above or equal to 25 Mq, this diagram is 
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sensitive to both mass loss and mixing. Since WN sub- 
types are sensitive to the optical thickness of the winds, 
while Xg is sensitive to both mixing and mass loss, one 
should try to define the domains of this diagram where 
" slash stars" , LBV, WN late and WN early are each lo- 
cated on average. Ideally this should also be attempted for 
stars in different galaxies. For masses below 25 Mq, this 
diagram indicates up to which stage mixing is proceeding 
for stars of different luminosities, which is also a useful 
indication. 

The evolution of the (C+0)/He ratio as a function 
of luminosity is shown in Fig. I20L which is the key dia- 
gram for WC stars as shown by Smith & Maeder 1)1991(1 . 
The evolution for a given mass is upwards, the luminos- 
ity decreases as the stellar masses decreases. At the same 
time, as a result of mass loss and mixing, the (C-|-0)/He 
ratio goes up, since there are more and more products 
of He-burning at the stellar surface. Rotation drastically 
modifies the tracks of WC stars in this diagram, as well as 
the lifetimes. The tracks of the models of 85 and 120 Mq 
are shifted to lower luminosities by rotation. As there is a 
tight mass-luminosity relation for He, C, O stars (Maeder 
I1983|l . this is the result of the smaller mass at the en- 
trance of the WC phase for models with rotation. The 
smaller mass is due mainly to the much longer WN phase 
of rotating models. (As a minor remark, we see that the 
(C-|-0)/He ratio goes slightly higher in models without ro- 
tation, this is due to the fact that the He- and C-burning 
proceeds a bit farther in non rotating models that have a 
higher mass). For the models of 60 Mq, there is a small 
difference at the beginning of the WC phase and this dif- 
ference quickly disappears. The 40 Mq models enter the 
WC phase only if rotation is significant. In a previous work 
(Smith & Maeder ri991|l . we have shown that metallicity Z 
strongly shifts the tracks in this diagram. Thus, we now see 
that rotation is also producing great shifts. The location 
of WC stars of different subtypes and composition should 
be attempted; it would be particularly worth considering 
WC stars of galaxies of different Z. 

8. Conclusion 

The good agreement obtained here between the observed 
WR star populations and the predictions of the rotating 
stellar models reinforces the general conclusion, already 
obtained in previous works, namely that rotation is a key 
physical ingredient of massive star evolution. Population 
synthesis for galaxies and starbursts should also be based 
on stellar models with rotation. The same is true for nu- 
cleosynthetic models. 

Rotation, in addition to making the models more real- 
istic, opens new ways to explore very interesting questions. 
Among them let we cite the possible relations between ro- 
tation and the ring nebulae observed around some WR 
stars, the relation between rotation and the LBV phe- 
nomenon, the rotation of pulsars and the link to 7-ray 
bursts. 
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